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In a recent paper Sanpera et al. have shown [1], that for the simplest binary 
composite systems any density matrix can be described in terms of only product 
vectors. The purpose of this note is to show that posibillity of decomposing any 
state as pseudomixtures does not depend on dimension of the subsystems. 

In a composite quantum system one can find states that involve different degrees of 
correlation between subsystems. From the physical point of view the most important 
characterization of such states is related to the nature of these correlations. Separable 
states involve correlations which can be explained by classical models [2]. The notion 
of inseparability or entanglement is related to specific relations existing only in pure 
quantum systems. Recently a general characterization of separable states was given in 
the language of vectors from natural cone of the Tomita-Takesaki theory, and quantum 
origin of this notion was explained [3] . 

Quantum character of entaglement plays a crucial role in quantum communication 
[4], cryptography [5], and quantum computation [6]. A state of a composite quantum 
system is inseparable if it cannot by represented as a convex mixture of tensor products 
of states of its subsystems (p ^ J2u)1'^P'^^ ® P^^^ otherwise the state is separable). 
This definition is very hard to handle. Therefore another characterizations are needed. 
Peres [7] and the Horodeckis [8] have obtained a characterization of separable states for 
systems described by Hilbert spaces with dimensions 2 x 2 or 2 x 3. Sanpera et al. 
have shown that for a binary system of dimensions 2x2 any state can be described as 
a linear combination of separable states called pseudomixtures. 

In this paper we show that any state of any composite quantum system described 
by Hilbert spaces with finite dimensions can be given in a form of a pseudomixture. In 
order to show this we use a simple modification of the definition of a separable state 
[9]: 

Fact 1 Any density matrix p is separable iff it is of the form: 

P=^PaPa (1) 

a 

where Pa > 0, J2a Pa — ^ Pa are projectors on simple-tensor vectors. 
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Let us consider two physical systems described by Hilbert spaces Tii and 7^2- A com- 
posite quantum system is described by the Hilbert space Hi (^0.2- We need some facts 
about self-adjoint Hilbert-Schmidt operators (p = p* and p € BusCHi ® H2 ))■ 

Proposition 1 For any self-adjoint Hilbert-Schmidt operator p and any simple-tensor 
basis {| Bj ig) fj)}ij there exists an orthogonal decomposition of p into a diagonal com- 
ponent A and an off-diagonal component H: 

p = A + H and {A,H)hs=0, (2) 

where {A, H)us = Ty AH is a scalar product in BhsCHi ® 7^2) • 

Proof: Let p = ^fPtPt bc' the spectral decoposition of p. Then we have: 

t t ij kl 

= X] Pt'^ij'^ki I e» ® fj){ek fi I 

ijklt 

= "^Ptlaijf \ei^ fj){ei | + ^ PtaljO^ki I ^ fj) (efe O /; | 

= A + H (3) 

and 
AH = 

«) /„) (e™ «) /„ I X PsaljUli I (g) /j-) (cfe (g. /j | 

m,n,t ijik,jjtl,s 
= Y Pt\ "mn 1^ P^^ij^i I ^™ ® /") ® /n 1 1 ei O /j> (efc <8) /i I 

m,n,t i^k,i^l,s 

PtP^ I I ® (^fi ® I • (4) 

Therefore Tr^iJ = 0. □ 

Proposition 2 If for any decomposition (2) all diagonal components A are equal to 0, 
then p is equal to 0. 

Proof: If we have a decomposition (2) associated with a basis {| ® fj)}ij, then any 
other decomposition can be obtained by a unitary transformations Ci — J2k''^ikek and 
fj = '^jlfl ill tlie following way: 



mn 



ij kl h mn k-i I2 

= I ^C^'^ljUil^i^jh)eki ® fh)(YC^'^mnUmk2Vnh)ek2 ^ fh I 
kih ij k2l2 

= X! X! '^ij'^mnUikiVjl^Umk2Vnl2 I hi <^ fli){ek2 ^ fh \ ■ (5) 

kilik2kil2 ijmn 
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It follows that the operators with the decomposition associated with the basis 
{| e/j (8) fi)}ki (which will be described by A and H) can be given as a function of a*^ 
(these coefficients will be referred to as the coefficients of spectral projectors in the basis 

-4 = ^ ^ aljal^^UikVjiUmkVni I efc ® fi){ek (8) /; | (6) 

t k,l ijmn 

and 

H = ^Pt X] m aljai,^Uik,Vji,Urak2Vni2 I 4i ^ fh){ek2 ^fhl- (7) 

t fci,ii,fc27^fci,^27^^i ij'rnn 

Putting Pijmn — ^tPi^ij^mn using Hcrmiticity of p we obtain: 

)\ek®fi){ek®fi\ (8) 

i>m,j>n 

and 

■) I 4i ®/h>(efe ^/;2 I • (9) 

/ci Ji,fc27^fei)^27^^i i'>m,j>n 

Assume that all operators A = for any decomposition (2). It means that 

■) = 0, (10) 

i>m,j>n 

where Ui is the set of all unitary transformations on Wj. We can choose unitary transfor- 
mations which in the fco-th column has only two nonzero elements Uaka = {^aio +5amo ) 
and analogously = '^{^bjo + ^buo)- Foi' such u and v we have: 

^SR(/3womonJ = 0. (11) 
Next, choosing Uako = -^{i^aia + ^ama) and V as before, we have 

^5?(iA„j„™„„J = ^S(/3z„j„m„nJ = 0. (12) 

Then Pi^j^rn„no = for all io rrio and jo ^ no, so p = 0. □ 

Proposition 3 If p is a density matrix then a norm of the operator H is less than ^ 

and its trace is equal to 0. 

Proof: We can write the operator H in the following form: 

^= E Pt'^lAi\<^i^fj){^k^fi\= E P^^m- (13) 

ijtk,jjtl,t i>k,j>l,t 
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The operator G^jj^i has two eigenvalues | al^a^i \ and — | a\ja\.i \. Since 

< (I <■ 1-14/ I)' =1 I 4/1' - 2 I I, (14) 

then 

l<"Ll<^(l<l'+l"Ll')- (15) 
Using this we can estimate the norm of the operator H: 

\\H\\ = II E PtGUi\\< E p*iiGUii= E f^*K-4;i 

i>k,j>l,t i>k,j>l,t i>k,j>l,t 

i>k,j>l,t t ijkl 

The last equality results from J2tPt ~ 1 Sij I ^ijl' ~ 1- Tr^ce of this operator 
equals 0: 

1 = Trp=TriJ2Pt\<^i/\ei®fj){ei^fj\)+TrH 

= Y^Pt I 41' Tr (I ei ® /,)(e, ®f,\)+TvH^ J^P^ I + Tr 

= E^'*EKi' + ^^ = i + '^^- (17) 

□ 

Let /9 be any density matrix from Bhs('Hi ®'H2)- For this matrix we can choose such 

a decomposition (2) that Tr has a maximal value. This maximum exists because we 
take supremum over U\®U2, so over a compact set, of a continuous function of u and 

iti2 = E(Ep*iE"l-^^'^^^'i')', (18) 

k,l t ij 

Therefore there exists a basis realising this maximum, namely {| <S) fj)}ij' 

P = E Ptol/3*lhe?®/°)(eO®/0|+ E PtoPlj^i\e1®f'^){el<^ff\ 

i—k,j—l,to i^kjj^l^to 

= A{1) + H{1) and TrA(l)2 = sup Tr A^, (19) 

where A describe the set of all diagonal components in simple-tensor decomposition, 
Pt^ and arc the eigenvalues of p and the coefficients of spectral projectors in a basis 
{| e° ® fj )}ij, respectively. 

Prom Proposition 3 we know that ||iJ(l)|| < ^. The operator ^4(1) is a linear (in 
the first step even convex) combination of simple tensors, so we can think about this 
decomposition as separating out the simple-tensor part. Since H{1) is self-adjoint (as 
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a sum of self-adjoint operators Glji~i) we can continue the separation procedure. We 
again choose a decomposition with maximum Tr which is associated with the basis 

+ E ''swPS'^pS^ I 4 ® fj){el ® fl 1= ^(2) + H{2) (20) 

i#fe,i#;,s(i) 

After n steps we obtain: 

n 

p = Y,M^ + H{n) (21) 

i=l 

Since the decomposition is orthogonal, we have 

TrH'^{n-l) = TvA'^{n)+ TvH'^in) (22) 

and 

^ I o,(„_i)|^ = E(E««W I + E I ''^Wi' > E I (23) 

s(n— 1) i,j s(n) s(n) s(n) 

The equality Tr H^{n — 1) = Tr i7^(n) can be obtained iff A{n) = and, therefore, by 
Proposition 2 only for H(n — I) = 0. Then a sum of squares of eigenvalues of H(n) is 
less then a sum of squares of eigenvalues of H{n — 1). This means that lim„ ||ff(n)|| = 
0. In the limit of these procedures we obtain a mixture of simple-tensor projectors 
{I eg)/") (e"(g)/" \}iju, not necessarily orthogonal to one another, with real coefficients. 
If we collect positive and negative coefficients separetely we obtain a pseudomixture: 

p = ap^«f(+) -fop^'^P^-), (24) 

where a,b > 0, a = 1 + b, and and arc separable states. For the simplest 

binary composite systems, it was shown [1] that by using geometrical properties typical 
of this dimension it is possible to determine cardinality of and i.e. the 

smallest a in the decomposition (1). Our procedure does not give such an information. 
The decomposition (24), similarly to other decompositions of this kind, is not unique. 
The most interesting pseudomixture is the one for which b is minimal. Our procedure 
does not determine such a pseudomixture. 

Since in our procedure we use only the properties of the Hilbert- Schmidt scalar 
product, it seems that it can be transferred also to composite systems described by 
infitc-dimcnsional Hilbert spaces. The language of vectors from natural cone of the 
Tomita-Takesaki theory is the most natural one for such a problem. The infinite- 
dimensional case will be discussed elsewhere [10] . 

I would like to thank for the hospitality to Centrum Leo Apostel in Brussels where 
a part of this work was done. I am indebted to W.A.Majewski and M.Czachor for 
comments. The work was supported by the Polish-Flemish Grant No. 007. 
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